Известия Томского политехнического университета. 2004. Т. 307. № 2 


УДК 622.831.232 

МОДЕЛИРОВАНИЕ ЦИКЛИЧЕСКОГО ХАРАКТЕРА ОБРУШЕНИЯ ГОРНЫХ ПОРОД ПРИ 
ПРОВЕДЕНИИ ВЫРАБОТКИ С ПОСЛЕДОВАТЕЛЬНЫМ НАКОПЛЕНИЕМ ПОВРЕЖДЕНИЙ 

Л.Д. Павлова, В.Н. Фрянов 

ГОУ ВПО "Сибирский государственный индустриальный университет", г. Новокузнецк 
Е-таіІ: гес1ог@$іЬ$іи.ги, кзОІІт.пкі.ги 

На основе метода конечных элементов разработана пространственная модель деформирования и разрушения горных пород для 
прогноза изменения напряженно-деформированного состояния углепородного массива и его прочностных свойств в процессе 
отработки угольного пласта. Механизм обрушения пород кровли моделируется в процессе пошагового движения очистного за- 
боя, с учетом влияния функции времени, анализа напряженного состояния в массиве горных пород, последовательного накоп- 
ления повреждений в породах для конечных элементов и обрушения тех из них, в которых накопленные повреждения превыси- 
ли критическое значение. 


Интенсификация процессов выемки угля в 
очистных забоях существенно влияет на геомеха- 
ническую ситуацию, связанную с напряженно-де- 
формированным состоянием (НДС) пород в окре- 
стности забоев. При этом важную роль в обоснова- 
нии параметров управления горным давлением 
имеет механизм формирования предельного нап- 
ряженного состояния в области интенсивного раз- 
рушения пород при сдвижении кровли, а также 
создание математических моделей, позволяющих 
прогнозировать изменение состояния углепород- 
ного массива при проведении горной выработки. 

Поскольку горные породы в процессе их дефор- 
мирования проявляют упругие, пластические и 
хрупкие свойства практически одновременно, то, ес- 
тественно, нельзя построить их математическое опи- 
сание на базе только одной из известных моделей. 

Для прогноза проявлений горного давления на 
основе оценки НДС массива горных пород метода- 
ми механики сплошных сред математическая мо- 
дель массива описывается тремя известными груп- 
пами уравнений: статическими уравнениями рав- 
новесия, кинематическими уравнениями совмест- 
ности и непрерывности деформаций и физически- 
ми уравнениями, связывающими компоненты тен- 
зоров напряжений и деформаций [1]. 

Решение уравнений осуществляется методом 
конечных элементов, который сводит задачу к ре- 
шению системы алгебраических уравнений путем 
минимизации функционала полной потенциаль- 
ной энергии среды. 

Кроме того, для методов, учитывающих воз- 
никновение вблизи выработки зоны разрушения, 
требуются реологические и прочностные характе- 
ристики горных пород, включая характеристики, 
определяющие закономерности деформирования 
пород в допредельном и запредельном состоянии. 

Разработка реологической модели для учета из- 
менения свойств угольного массива во времени. Для 
учета изменения свойств угольного массива во вре- 
мени используется реологическая модель, которая 
отражает свойство ползучести горных пород и ха- 
рактеризуется определенной зависимостью дефор- 
мации в какой-либо момент времени от всей пре- 
дыдущей истории нагружения. 


Согласно теории Ж. С. Ержанова, для описания 
поведения различных видов пород, свойства кото- 
рых во времени не изменяются, можно использо- 
вать разностное ядро типа Абеля. Если в момент 
времени 1 к телу приложено напряжение сг( г), 
действующее в течение времени ёт, то при исполь- 
зовании ядра ползучести типа Абеля [2] 

Ці-т) = 8(1- т)““, 


закон деформирования образца горной породы 
при длительно действующей постоянной нагрузке 
имеет вид [2] 


е(0 = 



8і 1 ~ а 


1 -а 


где е({) - деформация в момент времени {; Е - мо- 
дуль деформации; 5 и а(0<а<1) - характеристики 
ползучести, получаемые экспериментально. 


Параметр ползучести, характеризующий реоло- 
гические свойства угля, пород кровли и почвы, оп- 
ределяется в зависимости от литологического типа 
породы по корреляционным формулам [3]: 


Го, 102 -0,17 -10“ 5 Д 


<5 = 


0,20 -0,42 -10- 5 Е, 
3,38 -0,88 -10” 5 Е, 


0, 14 -2, МО' 5 Е, 


для песчанника ; 
для алевролита , 
для аргиллита , 
для угля, 


где Е - модуль деформации породы. 

Коэффициент а, учитывающий трещинова- 
тость угольного массива, также определяется экспе- 
риментально в зависимости от коэффициента кре- 
пости и расстояния между трещинами (а=0,4...0,8). 

В задачах механики подземных сооружений, в 
которых граничные условия и объемные силы при- 
нимаются не зависящими от времени, Б.З. Амусин 
и А.М. Линьков предложили использовать метод 
переменных модулей [4] . 

Согласно этому методу влияние времени учи- 
тывается путем замены деформационных характе- 
ристик массива временными функциями. В част- 
ности, модуль деформации пород Е можно предс- 
тавить как некоторую функцию времени Е, [4] 

Е = Е 
' 1 + *>(*)’ 
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где ср(1 ) - функция ползучести горных пород, учи- 
тывающая изменение свойств угольного массива 
во времени 


<р(і) = 


8і'~ а 
1 -а ' 


Временные функции для коэффициента Пуас- 
сона и модуля сдвига имеют вид [4]: 


ѵ,= 0,5- 


0,5-ѵ 

1-<р(0 


г 2(1 +ѵ) 

\ + 3<р(і)' 


Следует отметить, что метод переменных моду- 
лей предполагает замену деформационных харак- 
теристик временными функциями для всего уголь- 
ного массива, однако изменение свойств во време- 
ни необходимо учитывать только в зоне влияния 
горной выработки. 

Поэтому в данной работе для того, чтобы учи- 
тывать функции времени только в зоне влияния 
выработки вычисляется коэффициент концентра- 
ции упругих напряжений по формуле 



где т ос , - октаэдрические касательные напряжения в 
массиве, вмещающем горную выработку, вычис- 
ленные через компоненты нормальных напряже- 
ний { СГ СТ СТ- Г, Т -- Г- д } [1], 

т ос = 3 Ѵ (<т ѵ -°у) 2 +(оу -<х ) 2 +(<т, -<у х ) 2 +6(г^ +У +т2 ), (1) 

тѴ< - октаэдрические касательные напряжения в 
нетронутом массиве, вычисленные через компо- 
ненты нормальных напряжений {<т“<т^сг“г^г^г^} по 
формуле (1). 

Тогда функция ползучести с учетом коэффици- 
ента концентрации напряжений в зоне влияния гор- 
ной выработки определяется следующим образом 

<Р Ь ,=Ѵ,К- ( 2 ) 

Использование формулы (2) позволяет изменять 
модуль деформации ( Е) угля и пород в каждом ко- 
нечном элементе в зависимости от степени влияния 
на него выработанного пространства по формуле 



Таким образом, изменение свойств угля и пород 
позволяет учесть влияние функции времени при 
вычислении нелинейных параметров НДС массива 
горных пород и определении границ зоны блочно- 
го обрушения горных пород. 

Выбор критерия прочности для прогноза характе- 
ра разрушения горных пород. Критерии разрушения 
не следуют из уравнений равновесия и движения, 
они являются дополнительными условиями при 
решении задач о предельном равновесии. 

Прогноз несущей способности материала мож- 
но выполнить с помощью критериев локального 
разрушения или макроразрушения. Однако в ре- 
альных материалах и горных породах вследствие 
микропластических деформаций локальные пики 


напряжений перераспределяются и не влекут за со- 
бой разрушения тела. При разрушении в полях 
сжимающих напряжений даже бурное развитие од- 
ной или нескольких трещин не ведет к катастрофи- 
ческой потере несущей способности. В таких усло- 
виях более работоспособными оказываются крите- 
рии макроразрушения, описывающие потерю не- 
сущей способности некоторого объема материала. 

При моделировании механизма разрушения уг- 
ля и пород в окрестности горной выработки для 
определения предельных напряжений применяет- 
ся обобщенный критерий Кулона- Мизеса [5], в ко- 
тором условие разрушения в пространстве главных 
напряжений интерпретируется поверхностью вра- 
щения в виде усеченного конуса: 

к,=(Сс ІМ > + оф-^); 


где 

Ѵ2Сс^ ; СсЩ<р 

а + СсІ§ <р у/2 (сг с + СсІ§ ср) 


су от Уосі ~ нормальное и касательное октаэдрические 
напряжения; а р - прочность на растяжение (ст р <0); 
С - сцепление; (р - угол внутреннего трения. 

На основе анализа пространственного НДС с 
использованием критерия Кулона-Мизеса опреде- 
ляется коэффициент остаточной прочности угля 
или породы (ку 1) в каждом конечном элементе: 


к Р = 


_ ргесі 


ргесі 


сг, < <т„ 


< 7 „„, > СГ 


( 5 ) 


где г "У и с т" р а д - предельные октаэдрические напря- 
жения, вычисленные по критерию Кулона-Мизеса. 

Границы блоков обрушения подработанных по- 
род кровли выделяются по коэффициенту остаточ- 
ной прочности: значения ку 1 определяют зону 
блочного обрушения угля и пород. 

Моделирование механизма обрушения горных по- 
род с последовательным накоплением повреждений. 
Экспериментальными методами исследования ус- 
тановлено, что в зависимости от горно-геологичес- 
ких условий смещение пород кровли в окрестности 
горных выработок происходит в виде расслоения и 
последующего послойного сдвижения, ступенчато- 
го опускания блоков по пересекающим слои тре- 
щинам, а также плавного опускания без видимых 
разрывов сплошности. При этом система релакси- 
рует в два основных этапа: накопления микроско- 
пических повреждений и магистрального макрос- 
копического разрушения [6]. 

Второй этап наступает при достижении некото- 
рого интегрального предела микроразрушений, за- 
висящего от параметров НДС и свойств горных по- 
род. Поэтому для моделирования процесса обруше- 
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ния пород кровли при движении очистного забоя 
необходимо использовать механизм последователь- 
ного накопления микроскопических повреждений 
и определения критического значения предела де- 
зинтеграции горных пород, позволяющего выде- 
лить границу зоны обрушившихся пород. 

Механизм обрушения моделируется в процессе 
пошагового движения очистного забоя, с учетом 
влияния функции времени и изменения модуля де- 
формации угля и пород в каждом конечном эле- 
менте в зависимости от степени влияния на него 
выработанного пространства по формуле (3), ана- 
лиза изменений НДС в массиве горных пород по 
формулам (4, 5), последовательного накопления 
повреждений в породах для конечных элементов и 


обрушения тех из них, в которых накопленные 
повреждения превысили критическое значение. 

Процесс обрушения пород кровли носит цик- 
лический характер: можно выделить первичный 
цикл и циклы вторичных осадок. При ограничен- 
ных размерах области исследования цикличность 
процесса обрушения имитируется пошаговым 
продвижением очистного забоя до первичного об- 
рушения пород основной кровли и последующим 
обрушением при вторичных осадках кровли посре- 
дством перехода к новой системе координат путем 
параллельного сдвига по оси ОХ на границу разру- 
шенных пород при первичном обрушении. 

На основе разработанного алгоритма проведено 
прогнозирование параметров разрушения пород 
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Рис. 1. Изолинии коэффициента остаточной прочности пород при движении очистного забоя от монтажной камеры на: а) 10 м; 
б) 30 м; в) 50 м 
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Расстояние от монтажной камеры, м 

Рис. 2. Эпюры опорного давления для угольного пласта с учетом функции времени и движения очистного забоя 


непосредственной кровли в процессе пошагового 
движения забоя по пласту 26 а при отработке свиты 
угольных пластов в условиях ОАО "Шахта "Еса- 
ульская". Порядок отработки пластов нисходящий, 
верхний пласт 29 а отработан, поэтому породы вы- 
ше отработанного пласта рассматриваются как сы- 
пучий материал в виде пригрузки. Свойства пород 
между отработанным верхним пластом 29 а и отра- 
батываемым нижним пластом 26 а приняты соглас- 
но геологическим отчетам. 

Как показано на рис. 1, в виде изолиний коэф- 
фициента остаточной прочности пород кровли, 
блоки пород, расположенные непосредственно над 
выработанным пространством в слоях подработан- 
ной толщи, обрушаются, образуя свод уступной 
формы. Значение коэффициента остаточной проч- 
ности равное единице (к= 1) определяет границу 
зоны обрушенных пород. Суммарная мощность 
обрушившихся слоев пород достигает Юм при от- 
работке угольного пласта 26 а мощностью 2 м. Рас- 
четная высота зоны обрушения соответствует фак- 
тической высоте, равной от 3-х до 5-ти мощностей 
вынимаемого пласта. Угол обрушения пород в вы- 
работанном пространстве составляет 50.. .90°. По 
мере подвигания забоя лавы в движение приходят 
новые участки подработанной толщи. Увеличение 
площади выработанного пространства приводит к 
изменению НДС и пластических свойств пород, 


что способствует увеличению области сдвижения и 
формированию зоны дезинтеграции (рис. 1, а, б, в). 

Следует отметить, что характер развития про- 
цесса разрушения оказывает существенное влия- 
ние на распределение напряжений вблизи вырабо- 
ток, на параметры опорного горного давления, не- 
упругих деформаций и величины смещений, а так- 
же на пучение и разрушение пород почвы и кровли. 
До первичного обрушения пород непосредствен- 
ной кровли на границе угольного массива всегда 
четко выражена область повышенных напряжений 
(рис. 2). После первичного обрушения пород кров- 
ли напряжения в угольном массиве снижаются на 
30.. .40 %. 

Полученные результаты моделирования по ха- 
рактеру процесса разрушения горных пород и ве- 
личине максимальных напряжений качественно 
соответствуют данным натурных наблюдений по 
распределению напряжений в угольном массиве, 
приведенным в работах [7, 8]. 

Таким образом, разработанная пространствен- 
ная модель деформирования и разрушения горных 
пород позволяет на основе расчета методом конеч- 
ных элементов параметров НДС углепородного 
массива, использования обобщенного критерия 
прочности и временных функций адекватно опи- 
сывать механизм обрушения пород кровли при 
движении очистного забоя. 
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